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Using Ewald's summation method we investigate collective electronic excitations (excitons) of 
ultracold atoms in parallel planar optical lattices including long range interactions. The exciton 
dispersion relation can then be suitably rewritten and efSciently calculated for long range resonance 
dipole-dipole interactions. Such in-plane excitons resonantly couple for two identical optical lattices, 
with an energy transfer strength decreasing exponentially with the distance between the lattices. 
This allows a restriction of the transfer to neighboring planes and gives rise to excitons delocalized 
between the lattices. In general equivalent results will hold for any planar system containing lattice 
layers of optically active and dipolar materials. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices can be considered 
as one of the most attractive quantum systems that oc- 
cupied a wide range of theoretical and experimental re- 
searches [ij. The controlability of the different system 
parameters makes it a test set-up for many condensed 
matter effects and models, and opened the door toward 
answering fundamental questions and getting a deep un- 
derstanding of condensed matter phenomena 01 • Most 
recent research is restricted to atoms in their lowest hy- 
perfine manifold. As a key phenomenon a quantum phase 
transition from a superfiuid to a Mott insulator phase can 
be achieved by controlling the depth of the optical lattice 
potential in changing the laser intensity [3, Q . 

Electronic excitations for ultracold atoms in optical 
lattices, which posses an energy many orders of magni- 
tude bigger than the kinetic temperature of the atoms, 
are generally avoided in this context. Due to their short 
life time and the momentum recoil in emission or ab- 
sorption processes of a photon, this energy can heat the 
system and destroy the stability of the Mott insulator 
phase. This view is limited by considering electronic ex- 
citations as independent entities that localized at each 
optical lattice site, which is the case for time-of-flight 
observation techniques [5] , and also for metrology in us- 
ing clock states 6]. These considerations, which rest on 
treating the optical lattice only as a tool to localized the 
independent atoms at a fixed zones in space, missed the 
novelty of the formation of collective electronic effects 
that can be induced through the direct electrostatic inter- 
actions among atoms at different sites and in exploiting 
the lattice symmetry. 

In previous work we generalized optical lattice dynam- 
ics to include collective electronic excitations of the atoms 
For sufficient density such an excitation gets de- 
localized in the lattice due to resonance dipole-dipole in- 
teractions, where energy transfer among atoms at differ- 
ent lattice sites and in momentum space is represented 
by a wave that propagates in the lattice with a fixed 
wave vector. It can be considered as a quasi-particle 



termed exciton [lO|, . The exciton quasi- momentum is 
distributed over the whole lattice atoms so that at the 
emission or absorption of a photon involving an exciton 
the whole lattice recoils and not a single atom. There- 
fore excitons do not necessarily destroy to the quantum 
phases and in particular the Mott phase as they exhibit 
a suppressed recoil effect. Interestingly the life time of 
excitons in low dimensional optical lattices as calculated 
recently can be very short due to superradiance for 
some polarizations and wave vectors or metastable for 
others. This fact provides a strong tool that allows the 
controlability of excitation life times in optical lattices. 

In the present paper we investigate electronic excita- 
tions for ultracold atoms in planar optical lattices in an 
improved exciton picture. While in previous works we 
neglected the long range part of the dipole-dipole inter- 
action in a nearest neighbor approximation [7|-|9[ , here we 
explicitly calculate the corrections from long range inter- 
actions. In addition we account for the anisotropic nature 
of excitons for a fixed but nonorthogonal direction of the 
excitation transition dipole. The explicit calculation can 
be performed using the Ewald's method [14J. 

Once the eigenstates for a single plane are found the 
treatment can be extended to the case of two identical 
parallel planes. Excitons in the two planes will couple 
and lead to resonant energy transfer. Even though the 
resonance dipole-dipole interaction among two atoms in 
each plane has an inverse cube dependence on the dis- 
tance, a summation over the whole plane leads to an ex- 
ponentially decrease of the effective exciton-exciton cou- 
pling. At small distances one can form hybrid excitons 
delocalized between the two lattices planes. The results 
can be generalized to a finite number of identical planar 
optical lattices to describe a quasi 3D setup. 

The paper is organized as follows: in section 2 we 
present hybrid excitons for a system of two parallel and 
identical planar optical lattices. The exciton dispersions 
are calculated in sections 3 in using the long range reso- 
nance dipole-dipole interaction and in applying Ewald's 
method. This section includes the calculation details of 
the exciton matrix elements for inter-lattices. The con- 
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elusions arc given in section 4. The appendix includes the 
calculation of the intra-lattice exciton matrix elements. 



to get the full free exciton Hamiltonian 



(6) 



II. EXCITONS IN PLANAR OPTICAL 
LATTICES 

We consider two parallel identical planar optical lat- 
tices at distance b and lattice constant a filled with atoms 
in the Mott insulator phase with one atom per site as in 
figure (1). The atoms are two-level systems with on-site 
electronic transition energy Ea- An on-site electronic ex- 
citation couples to other sites via resonant dipole-dipole 
interaction, where an excitation at one site decays and 
other site gets excited. The energy transfer can be among 
sites in the same plane or among different planes. The 
electronic excitation Hamiltonian is given by 



(1) 



The indexes (n, m) run over the position of the in-plane 
sites, and (a, /3) over the two planes which are denoted 
by (1,2). Here Bl^^, Bna are the creation and annihi- 
lation operators of an electronic excitation at plane a 
and site n. They are taken to be boson operators with 
the commutation relation [Bna, bI^^] = SnmSa^- Using 
boson operators for two-level atoms is an approximation 
that hold for a single excitation or for low density of exci- 
tations, where we neglect saturation effects in which two 
excitations trying to excite the same site. The parameter 
"^nm is for the resonance dipole-dipole interaction among 
an excitation at site (n, a) and site (m, (3). 

Defined by the lattice symmetry, delocalized eigen- 
states can be represented as waves that propagate in the 
lattice with in-plane wave vector k and are denoted exci- 
tons. Hence the following transformation into the exciton 
basis formally diagonalizes the Hamiltonian 



Br. 



(2) 



where is the number of in-plane sites. The Hamilto- 
nian casts into 

H = ^Ea bI^B^^ + ^ J"^(k) Bl^Bup, (3) 



where we defined the exciton dynamical matrix 



R 



ik-R 



(4) 



The interaction parameter is a function of the distance 
between the two sites, that is = J"'^ (n — m) , where 
we defined J"''(n - m) = J"''(R), with R = n - m. 

As we have two identical lattices, the Hamiltonian can 
also easily be diagonalized in the lattice indexes in using 
the transformation 



-Bkzb = 



Ski ± B] 



k2 



(5) 



with {i/ = ±), and the exciton eigenenergies are 

(k) = Ea + J(k) ± J'(k), (7) 

where wc defined inside the same lattice the coupling 

parameter 

J(k) = J"(k) = J22(k), (8) 

and between different lattices the coupling parameter 

J'(k) = Ji2(k) = J2i(k). (9) 

Our main target now is to calculate the exciton dynami- 
cal matrices J(k) and J'(k). 
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FIG. 1: Two parallel planar optical lattices of lattice constant 
a, and separated by distance 6, with one atom per site. The in- 
plane wave vector, the in-plane and normal transition dipole 
components are plotted. 



III. EXCITON DISPERSIONS 

Lets us now calculate the exciton energy dispersion 
explicitly for dipole-dipole interaction between two sites 
separated by a distance R, which are given by 



J(R) = 



1 



IRplMp 



3(/i • R)2 



47reo 



|R|5 



(10) 



where /lx is the electronic excitation transition dipole, 
with the components /j, = {fix, f^y, fiz)- The general 
distance between two sites is R = {lxa,lya,lzb), where 
the lattice is taken to be in the {x — y) plane. We as- 
sume large square optical lattices with the number of 
in-plane sites N and area S = Na^. Hence, lx,ly = 
0, ±1, ±2, • • • , ±y/N /2, with the limit of iV » 1. As we 
consider here only two lattices, we have Iz = 0, 1, where 
one plane with z = and the other at z = b. The inter- 
action parameter is written as 



J(R) 



1 



47reo 



where = 

A,(R) = 



X, 2/, z), with 
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where Ix = alx, ly = aly and Iz = biz- In Fourier momen- 
tum space we get 



where 



R 



(13) 



(14) 



We have k = {kx,ky), which takes the values kx,y = 
^p, with p = 0, ±1, ±2, • • • , ±VN/2. 

From this point on we spht the calculations into two 
parts, one for in-plane interactions, and the other for in- 
teractions among the two lattices. Here we concentrate in 
the exciton dynamical matrix elements for inter-lattices, 
while the intra-lattice case appears in the appendix. 

For interactions between the two optical lattices we 
have = 1. We define the function 



5(k) = J2 



Then we have the diagonal elements 

;,2 



5/2 ■ 



Dyyik) - (2 

DzzO^) = 

and the off-diagonal elements 

Dxy{k) = D*{k) = 3 
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52 j S{k), 

- 2b^\ 5(k), (16) 
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dkxdk, 



-5(k), 



Dxzik) = Dl.ik) = Sib^Sik), 
Dyzik) = D:y{k)=m-^S{k). 



(17) 



We proceed in the calculation of the function S{k). 
Due to oscillations the series terms are slowly converging. 
Then we use Ewald's method to convert the oscillating 
series into a one with exponentially decaying terms [l3 |. 
and then we can concentrate in the first terms that rep- 
resent small wave number or long wave length excitons. 

In using the relation 



(18) 



we obtain 
S{k) = 



(19) 



We use the identity, for [j — x,y), 



'7rn+- 



to get 



S{k) 



3a2 



where 



r„„, = Tin 



dt ti/2e-^' 
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(20) 



(21) 



By applying the result 



/ dt ^fte-'^'e-^' 
Jo 



-ct^-h/t _ l + 2Vch /^^_2V^ 
2c 



^1 . (22) 



(23) 



— e 



(15) we get finally 



(24) 



The result of equation (|24|) is the main result here, as 
we got a new series with terms that decay exponentially, 
rather than the one we started with of equation (jl5l) with 
oscillating terms. Besides giving an intuitive picture this 
form is particularly useful for numerical calculations as 
the coefficients are real and decay fast. We can also easily 
treat some interesting limits in this form. 

As an example at this point on we concentrate on long 
wave length excitons with small k. In this limit we have 
fca <C 1, and hence we can keep only the first term of 
the series, the one with (n, m = 0), which can be easily 
justified numerically, then we get Fqo — ka/2, to obtain 



5(k) 



2tt 
3^ 



(1 + kb) e-'"'. 



(25) 



We calculate now the exciton dynamical matrix terms. 
We obtain the diagonal elements 

a'' k 

D (k) - —^e-^^ 



^yy 

DzzO^) = - 

and the off-diagonal elements 

Dxy{k) = D;,{k) 

Dxzik) = Dl^ik) 
Dyzik) = Dlyik) 



27r, 



-ke' 



-kb 



(26) 



-I — ttKxG 

a2 " 



-kb 



-kb 



(27) 
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where k = ^ + k'^. 
Finally we get 



J'(k) 



1 27r 
47reo o,^ 



/i^ \ ka e 



-kb 



(28) 



where we defined /x = (fipiiz) with = (ii^jfJ-y), and 
the unit vector k = k/fc. 

The in-plane calculations for the exciton dynamical 
matrix elements, which are given in the appendix, yields 



J(k) 



F 



(29) 



where calculating the first terms of the series (|A21I) gives 
F ~ 9/2. Note that the calculation for the limit ka <C 
1 in considering only the nearest neighbor interactions 
gives F = A. 

Here we present the results for typical optical lattice 
numbers, where the transition energy is Ea = 1 eV, the 
lattice constant is a = 1000 A, the transition dipole is 
fi = I eA, the interlattices distance is 6 = 10a, and the 
wave number is taken to be ka = 10^"^. We define Jq — 
^jf^^i , which is Jo ~ 1.44 x 10"^ eV. In figures (2- 
4) we plot the scaled transfer parameter, J'(k)/Jo, as 
a function of the angle between the in-plane wave vector 

and the transition dipole, (/), where • k^ = cos0, 

for different direction of the transition dipole, where fi — 
^{sm9,0,cos9). Figure (2. a) is for 9 = 0, figure (2.b) 
for 9 = 7r/6, figure (3. a) for 9 = tt/5, figure (3.b) for 
9 = 7r/4, figure (4. a) for 9 = tt/B, and figure (4.b) for 
6* = 7r/2. The plots show strong anisotropic effect. 
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FIG. 2: The scaled transfer parameter, J'(k)/Jo vs. the an- 
gle between the wave vector and the transition dipole, 0, for 
transition dipole direction of: (a) 6 = 0, and (b) 6 = n/6. 



IV. CONCLUSIONS 

Here we give a discussion of the results we obtained and 
present some conclusions. For the system of two identical 
and parallel optical lattices, an exciton that propagates in 
one lattice with a fixed wave vector k can hop to the other 
one due to dipole-dipole interactions. Such a coupling 
can form hybrid excitons among the two optical lattices, 
where we get two hybrid exciton modes, symmetric and 
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FIG. 3: The scaled transfer parameter, J'(k)/Jo vs. the an- 
gle between the wave vector and the transition dipole, </!>, for 
transition dipole direction of: (a) 9 = n/5, and (b) 9 = 7r/4. 
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FIG. 4: The scaled transfer parameter, J'(k)/Jo vs. the an- 
gle between the wave vector and the transition dipole, </!>, for 
transition dipole direction of: (a) 9 = n/3, and (b) 9 = -k 12. 



antisymmetric modes, which are separated by a splitting 
energy of A(k) = (k) - E-^[\^) = 2J'(k), where we 
found the result 



A(k) 



1 



ka e 



-kb 



(30) 



The result shows strong anisotropic effect in changing the 
direction of k relative to the in-plane transition dipole 
as seen in figures (2-4). Furthermore, polarization 
splitting appears at = between the normal component 
/iz and the in-plane component /xy, which equals ^^f^s ■ 

The most important feature here is the exponential de- 
cay behavior as a function of the inter-lattices distance 
of the exciton coupling among the two optical lattices. 
Even though the resonance dipole-dipole interaction has 
inverse cubic dependence of the distance between two lo- 
calized dipoles, the appearance of excitons induce expo- 
nentially decrease dependence for the exciton collective 
dipole among the two optical lattices. 

For the limit of & ^ a we can neglect the hopping of ex- 
citons and it is a good approximation to consider the two 
lattices independently. While for & ~ a we get the above 
hybrid symmetric and antisymmetric excitons. Note that 
the antisymmetric excitons are dark and metastable with 
long radiative life time, while the symmetric excitons are 
bright with damping rate two times larger than an ex- 
citon in a single optical lattice. Inside a cavity, in the 
strong coupling regime, the antisymmetric excitons de- 
couple to the cavity photons, while the symmetric ex- 
citons are coherently mixed with the cavity photons to 
form cavity polaritons. We treated similar case in our 
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previous works 043- 

The present set-up can be generalized to include finite 
number of parallel planar optical lattices. Then we get 
hybrid exciton modes as the number of the optical lat- 
tices, where part of the modes are dark and others bright. 
The above important result can be adopted in order to 
simplify the treatment by assuming hopping of excitons 
only among nearest neighbor optical lattices. Note that 
the hopping can be controlled by changing the distance 
between the optical lattices. The results hold for any 
similar system of two or multi layers of optically active 
lattices, e.g., for monolayers of molecular crystals, or lat- 
tices of quantum dotes. 

The results of the present paper can be adopted, with 
the appropriate modifications, for dipolar bosonic ultra- 
cold quantum gases in optical lattices, which interact 
via dipole-dipole interactions [H, [l^. One candidate 
for such a set-up is heteronuclear polar molecules that 
have permanent dipoles. The molecules need to be in 
their lowest rovibrational states, and an external electric 
field is applied to orient the molecules along a fixed di- 
rection. The other candidate is atomic species with a 
large magnetic moment that experience magnetic dipole- 
dipole- interactions, e.g., a BEG of ^^Cr, which have been 
realized recently in low dimensional optical lattices [l7j . 

The work was supported by the Austrian Science 
Funds (FWF), via the project (P21101). 



Appendix A: Intra-Lattice Exciton Dynamical 
Matrix 



Here we present the detail calculation for the exciton 
dynamical matrix elements for the in-plane interactions 
in using the long range resonance dipole-dipole interac- 
tions and in applying Ewald's method [14]. 

Now we take = in equation (fT2|) . We define 



5/2 



g»a(fe./.+fc,/,)^ (j = x,y), (Al) 



where the prime on the summation excludes the term 
with Ij = 0. The exciton diagonal dynamical matrix 
elements are 



Dyy(k) = ^2Sy{k)+S^{k), 
DM = S,ik)+Sy{k). 



(A2) 



111) yields 



(A3) 



Applying relation (f20l) yields 



Hoo +C30 



Sx{k) = ^E E COs(afc:r?a;) ( -^^j^T^ 



— 1 n— — oo 



(A4) 



We use now 



^ dte-'^'^e-''' =2^Ki{2^fch), (A5) 

where Kn{x) is the modified Bessel function of order n. 
We achieve the expression 



+ C30 +00 



A! 



S,{k) = E cos{akJ^)-j§^ 

l^—ln—-~oo ^ 

X {AtlK[{A':i)~K,{A':i)}, (A6) 
where the prime denotes derivative relative to l^, and 

Atl^ 21^ [nn+^y (A7) 

We use the relations 

2n 

Kn+i{x) = X„_i(x) + — X„(x), (A8) 

X 



to get 



+00 +00 



a: 



n— — 00 



/2 



X I^KoiAlD + K.iAll)]. (A9) 



We move to calculate the off-diagonal term. We define 

o a2 ' pia{k^l^+kyly) 



Using the relation ^TE\\ we get 

4 92 /.oo 



DM) 



' dkxdky 



In the following we concentrate in the calculation of 
Sxik), as the one for Syik) is similar. Using relation 



(All) 
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and relation pO)) gives 

+00 +00 



where T{n) is the Gamma function that takes r(2) = 1, 
to obtain the form 



Dxy{k) = 8 ^ ^ ( 7rn + — ^ j sin(a;2;A:^) 



ak. 



f 

Jo 



dt 



2 

afc„ \2 



(A12) 



^, ^. 8^1 Stt^ 4vr2 



(A17) 



and in using (|A5[) we obtain 



(ii). For = we get 



^^.(k)=4E E ^sin(a/.fc.)ir,(AX). (A13) ^^^^ = ^ yE^E^(-) 



n— — 00 



Here we going to calculate the result for the limit of 
long wavelength excitons, namely for fca <C 1. We start 
from the recent result of the off-diagonal term, in this 
limit we keep only the n — term, to get Dx.y{ka ^ 1) — >■ 
0. Note that for a; we have Kn{x +00) ^ \/~^^^^ 
and the series converges. 

We return now to the diagonal terms. In the limit 
of long wavelength excitons we consider the two main 
contributions: 

(i). For n = 0, in using cos(afc2:/a;) — > 1, we get 

Sx(n = 0) = o E \ ^o{akylx) + Ki{akylx) > , 
o aKylx J 

(A14) 

and using (|A8D yields 



+CX3 



Sx{n = 0) = o E kla^K2{akylx 



We use also the relation 



Kn{x < 1) 



r(n) (2 



2 \x 



(A15) 



(A16) 



X |i^o(27rn;,) + ^— -^i^i(27rnZ,)|4l8) 
where the prime exclude (tt, = 0). In using (jASp we have 



Sx{k = 0) = — — ^ n^K2{2TTnm). (A19) 



2 +°° 



Finally we get the result 

Sx{ka < 1) = Sy{ka < 1) ~ i^. 



where 



47r2 32^2 ^ 



(A20) 



(A21) 



For the diagonal dynamical matrix elements of excitons 
we obtain finally 



(A22) 
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